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1.  INTRODUCTION 

From  similarity  theory,  the  drag  and  thermal  ex¬ 
change  coefficients  (C[>,  CH,  and  Cr)  for  10m  above  the 
ocean  surface  depend  on  the  atmospheric  stability  pa¬ 
rameter  (z/L,,)  where  L„  the  Obukhov  length  scale  for  the 
atmospheric  surface  layer.  Cloud,  precipitation,  and 
ocean  mixed  layer  (OML)  thermodynamical  feedback 
affects  the  atmospheric  stability  parameter  in  two  ways, 
Clouds  reduce  the  incoming  solar  radiation  at  the  ocean 
surface  by  scattering  and  absorption,  which  cools  (rela¬ 
tively)  the  ocean  mixed  layer.  The  OML  cooling  lowers 
the  sensible  and  latent  heat  fluxes  from  the  ocean  sur¬ 
face,  which  increases  the  atmospheric  stability  parame¬ 
ter.  On  the  other  hand,  precipitation  from  the  clouds 
dilutes  the  surface  salinity,  stabilizing  the  upper  ocean, 
and  reducing  OML  deepening.  The  ocean  mixed  layer 
may  be  caused  to  shallow  if  the  downward  surface 
buoyancy  flux  is  sufficiently  enhanced  by  the  precipi¬ 
tation.  When  the  ocean  surface  receives  a  downward  net 
surface  heat  flux,  the  reduction  in  the  OML  depth  will 
increase  the  sea  surface  temperature  (SST)  by  concen¬ 
trating  the  net  radiation  plus  heat  fluxed  downward 
across  the  sea  surface  into  a  thinner  layer  (Fig. I).  The 
increase  of  SST  augments  the  sensible  and  latent  heat 
fluxes  from  the  ocean  surface,  which  decreases  the  at¬ 
mospheric  stability  parameter.  On  the  other  hand,  when 
the  ocean  surface  has  a  net.  surface  heat  loss,  the  re¬ 
duction  in  the  OML  depth  will  decrease  SST  due  to  the 
decrease  of  the  thermal  inertia.  The  decrease  of  SST  re¬ 
duces  the  sensible  and  latent  heat  fluxes  from  the  ocean 
surface,  which  increases  the  atmospheric  stability  pa¬ 


rameter.  The  variation  of  atmospheric  stability  parame¬ 
ter  in  turn  changes  the  drag  and  thermal  exchange 
coefficients.  In  this  study,  a  simple  coupled  eloud- 
prccipitation-ocean  mixed  layer  model  proposed  by  Chu 
ct  al.  (1990)  and  Chu  and  Garwood  (1991)  is  used  to 
compute  the  time  rate  change  of  drag  and  thermal  ex¬ 
change  coefficients  caused  by  this  mechanism. 


2.  DRAG  AND  THERMAL  EXCHANGE  COEF¬ 
FICIENTS 


From  the  Monin-Obukhov  similarity  theory,  the 
momentum,  heat,  and  moisture  fluxes  from  the  air-ocean 
interface  are  directly  related  to  the  mean  profiles  of 
wind,  temperature,  and  mixing  ratio: 


dz  L„ 


dTa 

dz 


T. 


^  4-  ,  ,  c  , 

dz  -  2  ^  La  } 


(1) 


where  <j>M,  <}>H,  <f>e  are  similarity  functions;  u.  is  the  at¬ 
mospheric  friction  velocity;  k  is  von  Karman  constant; 
T.  and  q.,  the  atmospheric  surface  layer  temperature  and 
humidity  scales,  represent  the  sensible  and  latent  heal 
fluxes  across  the  ocean  surface  (upward  positive), 

Qs=~  PacpaK“‘  T;  Qe  =  ~  Pa^u.q.  (2) 


where  pa,  c^,  L,  are  surface  air  density  (characteristic 
value),  specific  heat  of  air,  and  latent  heat  of 
vaporization.  For  a  near-neutral  atmospheric  surface 
layer  (small  z/LJ),  the  drag  coefficient  is  computed  by 


Cn  = 


[ln(f-)  +  ^]2 


(3) 


where  z,a  =  10m,  z0  is  the  roughness  length,  and 


d{zlLJ 


L/za-n-3-0 


from  the  Kansas  experiment  (Businger  et  al.,  1971).  The 
atmospheric  Obukhov  length  scale  Lt  is  defined  by 


La  =  - 


Kg(-#-+o.6ira-%-) 


(4) 
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where  7,  and  7,„  are  mean  air  and  virtual  temperatures. 

In  the  neutral  atmosphere,  z/L„  =  0,  the  drag  coeffi¬ 
cient  becomes 

C0.v  =  K2[!n(i^)]2  (5) 

If  bzl0/La  is  small  compared  to  ln(z,o/zo),  the  drag  coeffi¬ 
cient  for  a  non-neutral  atmosphere  is  related  to  the  drag 
coefficient  for  a  neutral  atmosphere  (C0.v)  by 

Co*Co*',1''^w)  {6) 

where  all  parameters  on  the  righthand-side  except  La  are 
independent  of  surface  heat  and  moisture  fluxes.  There¬ 
fore,  the  time  rate  of  change  of  drag  coefficient  due  to 
the  change  of  surface  heat  and  moisture  fluxes  is  influ¬ 
enced  only  by  the  change  of  L ^ 

I  SCD  ^  2bzl0IL‘  dLa 
CDN  dt  ln(r,o/aD)  dt 

Substitution  of  (4)  into  (7)  leads  to  a  expression  repres¬ 
enting  the  effect  of  change  of  surface  heat  and  moisture 
fluxes  on  the  time  rate  change  of  drag  coefficient: 


2^io*g  -  1 

ln(zi„/z0)7>.3  p“cP<- 


0.61  7- 


The  bulk  formulas  for  the  sensible  and  latent  heat  fluxes 


Pacpa^H^s  ~  2to)^  a>  Qe  ~  Pa^-v^-  *?10J ^'a 


well-mixed  turbulent  boundary  layer  which  exchanges 
heat  and  moisture  with  the  atmosphere  and  entrains 
water  from  deeper  layer.  The  heat  and  salinity  equations 
take  the  forms 


Kir=~ w‘(s ~ s-h) +{E~ Pr)s +As  {nb) 

where  is  the  mixed  layer  depth;  T_-„  and  S_h  are  tem¬ 
perature  and  salinity  at  the  base  of  the  entrainment  zone; 
E  is  the  surface  evaporation  rate;  Pr  is  the  precipitation 

rate;  Ar  and  As  are  horizontal  advection  for  temperature 
and  salinity,  respectively;  and 

F=Rb  -Rs  +  Qs+Qe  (13) 

is  the  net  heat  flux  at  the  air-ocean  interface,  where  R , 
and  R/,  are  solar  and  back  longwave  radiation  at  the 
ocean  surface.  The  entrainment  velocity  w,  is 
parameterized  by  Chu  and  Garwood  (1988)  as 

_ 'Ufi _  (14) 


p  =  (cyw.  -  c2Bhj 

is  the  surface  forcing  function;  C,  and  C2  are  tuning  co¬ 
efficients;  u*.  =  (pJpM),u Uj  is  the  water  surface  friction 
velocity;  and 

B  =  ^-  +  Pg(Pr-E)S  (15) 


where  C„  and  CE  are  the  10m  thermal  exchange  coeffi¬ 
cients.  If  the  mean  surface  wind  and  air  temperature  is 
determined  mostly  by  the  large  scale  atmospheric  mo¬ 
tion,  substitution  of  (9)  into  (8)  leads  to 

I  5C«  ST, 

-pr - =  M—*r  (1°) 

CDN  ot  dt 


where 


Here  we  assume  that 


JVflWrfln; 

Jd>w(0)rfln. 


1  dCD  i  dC„  1  dCF 
CDiV  dt  ~  CH  dt  ~  CE  dt  U  U 

If  we  use  7,as; Tlb—  7,— 288“A(,  and  Ui0^5mjs  the  M- 
values  is  roughly  estimated  as  A/«0.13/C-1, 

3.  OCEAN  MIXED  LAYER  DYNAMICS 

Most  models  that  include  thermodynamic  effects  in¬ 
clude  the  assumption  that  the  upper  ocean  layer  as  a 


is  the  downward  ocean  surface  buoyancy  flux.  The 
symbol  A.  is  a  Heaviside  function  of  P.  When  P  >  0, 
there  is  sufficient  turbulent  kinetic  energy  to  entrain  and 
mix  water  from  below  and  A„  =  1,  which  represents  the 
entrainment  regime.  The  time  rate  change  of  mixed  layer 
depth  is 


Here  w.*  is  the  mean  vertical  velocity  at  the  mixed  layer 
base.  When  P£  0,  there  is  insufficient  turbulent  kinetic 
energy  to  entrain  water  from  below,  and  w,  is  set  to  be 
zero,  i.e.,  A„  =  0.  This  is  called  ihe  shallowing  regime. 
The  mixed  layer  depth  is  calculated  from  7  =  0,  and  it 
is  directly  proportional  to  the  oceanic  Obukhov  length 
scale,  Lj. 

06*) 

Chu  and  Garwood  (1991)  discussed  the  effect  of 
cloud  amount  and  precipitation  perturbations  («',  Pr') 

on  the  OML  dynamics.  Three  parameters  were  found 
important: 
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PgSiPr^  E)  ^  sst  equation  for  the  shallowing  regime  is 


where  T  is  the  cloud-radiative  forcing  at  the  ocean  sur¬ 
face;  c  is  the  cloud-prccipitation  parameter,  representing 
precipitation  brought  on  by  unit  cloud  amount;  and  y  is 
the  surface  water  budget  index  representing  the  relative 
importance  of  surface  freshwater  influx  in  the  surface 
buoyancy  flux. 

Let  us  consider  the  effect  of  cloud  and  precipitation 
perturbation  on  SST  for  the  two  different  OML  regimes. 
After  mathematical  manipulation,  the  SST  perturbation 
caused  by  cloudiness  disturbance  for  the  entrainment 
regime  is 

4-(4r*’r>r,- 


Substitution  of  (19a,b)  into  (10)  leads  to 


for  the  entrainment  regime,  and 


for  the  shallowing  regime. 

In  determining  the  change  of  the  drag  coefficient, 
two  parameters  (,  y  are  crucial  when  OML  is  in  the 
shallowing  regime;  however,  <*,  T  are  important  when 
OML  is  in  entrainment  regime.  In  the  following,  we  only 
present  the  results  for  the  shallowing  regime. 


and  for  the  shallowing  regime  is 


where 


C  P^pJ^AT  18  day 

for  FslOOW'Vr2,  /i„~50m,  and  AT^3°K.  Here,  rr  is  the 
residence  time  for  SST  perturbation,  and  tc  is  the  time 
needed  for  OML  warming  up  to  the  deeper  water  tem¬ 
perature  by  F. 


4.  CHANGE  OF  DRAG  COEFFICIENT  DUE  TO 
HEAVY  PRECIPITATION 

Fig.2  shows  the  feedback  loop  among  cloud,  precipi¬ 
tation,  OML,  and  the  drag  and  thermal  exchange  coeffi¬ 
cients  ( CD ,  CE,  C„).  Cloud  and  precipitation  changes  the 
surface  buoyancy  flux,  which  in  turn  changes  SST 
through  the  variation  of  both  surface  net  heat  flux  and 
the  OML  depth.  The  SST  variation  changes  the  surface 
sensible  and  latent  heat  fluxes,  which  causes  the  modifi¬ 
cation  of  the  drag  and  thermal  exchange  coefficients. 

For  a  very  heavy  transient  precipitation,  such  as 
Pr  >  0.1  m/day. 


Fig.2.  Feedback  loop  among  cloud,  precipitation,  OML, 
and  the  drag  and  thermal  exchange  coefficients. 


Fig.3  shows  the  dependence  of  I ICDNdCnldt  on  the 
cloudiness  disturbance  n'  for  different  values  of  y  (0.1, 
0.5,  0.9,  2.5,  5.0)  at  i  =  0.5  m/day.  For  a  heavy  precipi¬ 
tation  case,  i.e.,  Pr-E> 0,  we  have  B>0  .  From  (17), 
y  >  1,  represents  the  upward  net  heat  flux  at  the  ocean 
surface  (ocean  losing  heat  at  the  surface);  and  y  <  1 ,  re¬ 
presents  the  downward  net  heat  flux  at  the  ocean  surface 
(ocean  gaining  heat  at  the  surface),  as  shown  in  Fig.4. 
Two  distinct  effects  of  cloudiness  disturbance  on  the  time 
rate  change  of  drag  coefficient  were  found  depending  on 
different  values  of  y.  For  y  <  I  (ocean  surface  losing 
heat),  the  atmospheric  surface  layer  is  destabilized  and 
the  drag  coefficient  is  increased;  however,  for  y  >  1 
(ocean  surface  gaining  heat),  the  atmospheric  surface 
layer  is  stabilized  and  the  drag  coefficient  is  decreased. 
For  y  =  5.0,  and  n'=0.5„  i/CD}idCDISt^  —  0.36/day, 
which  means  that  the  drag  coefficient  CD  will  be  reduced 
by  36%  during  l  day. 

TD/POL 


when  the  time  scale  for  the  precipitation  is  much  shorter 
than  rr  ,  the  SST  equation  for  the  entrainment  regime 
becomes 


I 


I 


Fig.3.  Dependence  of  I IC^CpJSt  on  the  cloudiness  dis¬ 
turbance  n'  for  different  values  of  y  (  0.1.  0.5,  0.9,  2.5, 
5.0). 


The  dependence  of  the  time  rate  change  of  the  drag 
coefTkicnt  on  n’  for  five  different  values  of  Z  is  shown  in 
Fig.4  for  y  ■*  0.1  and  in  Fig.5  for  y  ■>  5.0.  The 
cloud-precipitation-OML  effect  on  the  drag  coefficient 
strengthens  as  Z  increases.  If  the  precipitation  rate  from 
a  severe  storm  is  0.125  m/day,  and  during  the  storm  the 
cloudiness  only  changes  25%,  the  cloud-precipitation • 
coupling  coefficient  tscQ.Sm/day,  For  y  =  0, 1  (ocean 
surface  losing  heat)  and  Z*0.5mlday, 
llCDydCDldl*0.Q[  -  QAjday,  which  means  that  the  drag 
coefficient  increases  1-10%  during  I  day  (Fig.4).  For 
y  -  5.0  (ocean  surface  gaining  heat)  and  {ssO .5mlday. 

1  ICt>ndCDldt*  -0.2 /day,  for  n'i0.25,  which  means  that 
the  drag  coefficient  decreases  around  20%  during  I  day 
(Fig.5). 


Fig.4.  Dependence  of  1 /Cj><vdCp/dr  on  the  cloudiness  dis¬ 
turbance  ri  for  y-0.1  (ocean  losing  heat),  and  for:  (a) 
{  -  0.01  m/day,  (b)  Z  =  O.OSmjday  ,  (c)  {  «  0.1  m/day,  (d) 
0.25 m/day,  and  (c)  0.5 m/day. 


Fig.5.  Dependence  of  1  IC^dCpJtt  on  the  cloudiness  dis¬ 
turbance  rf  for  7  =  5.0  (ocean  gaining  heat),  and  for  (a) 
(  -  0.01m/<iry,  (b)  {  =  0.05 m/day  .  (c)  {  -  0.lm/<4rv.  (d) 
0.25 m/day,  and  (e)  0.5 m/Jay. 

5.  DISCUSSION 


(a)  In  this  study,  a  new  thermodynamic  feedback 
mechanism  among  cloud,  precipitation,  and  OML  is 
found  important  for  determining  the  drag  coefficient. 
This  process  should  be  considered  supplemental  to  the 
usually  dominating  factors  of  surface  roughness  on  the 
neutral  drag  coefficient. 

(b)  Heavy  precipitation  affects  the  drag  coefficient 
greatest,  especially  for  a  warm  ocean  surface  (i.e.,  ocean 
losing  heat  the  atmosphere).  The  drag  coefficient  can  in¬ 
crease  20-30%  within  I  day. 

(c)  It  is  not  the  purpose  of  this  first  study  to 
parameterize  the  drag  coefficient  In  detail.  In  the  future, 
we  will  use  more  sophisticated  large  eddy  simulation 
(LES)  model  for  both  atmosphere  and  ocean  to  do  fur¬ 
ther  research. 

(d)  Similar  to  the  drag  coefficient,  the  thermal  ex¬ 
change  coeffcients  Cu  and  C,  are  also  affected  by  this 
cloud-precipitation-  OML  feedback  mechanism. 
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